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Abstract 

We introduce the strange quark mass to the Sakai-Sugimoto model of holographic 
QCD. We compute mass shifts in the spectra of three-flavor baryons at the leading 
order in perturbation in quark masses. Comparison with experimental data shows 
an agreement only qualitatively. 



Oct 2009 



*e-mail: koji@riken.jp 
^e-mail: norihiro . iizuka@cern. ch 
^e-mail: ishiitk@riken.jp 
^e-mail: kadoh@riken.jp 



Contents 



3 



2 



1 



Introduction and summary 



A review of baryon mass shift 



Three flavors 



2 

8 



4 



Discussions 



12 



1 Introduction and summary 

Holographic QCD, a string theory realization of QCD in the large N c limit, offers a 
new approach to QCD. The Sakai-Sugimoto model [U E] is the most successful model 
of holographic QCD. Although the original proposal is a holographic dual of massless 
QCD, quark masses can be introduced by worldsheet instantons [3111] (see also [5j EJ [7] jE 
In the case of two flavors, mass shifts in hadron spectra due to the quark masses have 
been computed in [13], where the mass shift of nucleon was found to be consistent with 
lattice QCD results extrapolated to the chiral limit [HI [15j [16j [17]. In the situation that 
quarks are massless, baryons in three-flavor Sakai-Sugimoto model have been studied in 
|18j . However, in real QCD, where we have three flavors below chiral symmetry breaking 
scale, mass of the strange quark is heavy compared with up and down quarks. Hence it 
is important to deal with the third flavor as a massive quark, especially to compare with 
the experimental data of hadron spectra. 

In this short paper, we introduce the strange quark mass to the Sakai-Sugimoto model. 
We compute mass shifts in the spectra of three-flavor baryons, at the leading order in the 
expansion in quark masses. Our main result is shown in the formula (I3.13p . and the 
coefficients appearing there are listed in Table |2j where we used two input parameters m p 
and f n from experimental data to fix the parameters of the model, then comparison with 
experimental data can be meaningful. 

Baryons are identified with an instanton-like soliton of the Sakai-Sugimoto model pp. 
This soliton solution has been studied in detail in [19] (see also [20ll2T|[22|l23|[24l[25|l26]). 
The mass of baryons can be evaluated as a static energy of this nontrivial soliton of the 
gauge theory [TH [19] , and also a mass shift of the baryons due to the quark masses can 
be studied as an effect of the additional quark mass term in the action [13]. Such a kind 
of mass shift of baryons has been studied also in the context of chiral solitons (for three- 
flavor Skyrmions, see e.g. [271 [28], while for a review including other chiral soliton models, 
see e.g. [29]). However, in the Sakai-Sugimoto model, all the coefficients in the action are 
fixed by only two input parameters from the meson sector. This is a remarkable difference 
from many other chiral soliton mo delfl 

1 Another way to introduce the quark masses is a tachyon condensation in the D8/D8-branes [El [9j 



2 The existence of unnecessary modes whose counterparts are absent in QCD, is a long-standing 



HQlELlJQl. 
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In Section [21 we briefly review computation of the baryon mass shift in the case of 
two flavors [13], with giving a brief summary of necessary ingredients in this paper. We 
begin with the action of the Sakai-Sugimoto model, introduce the quark masses, give 
the baryon configuration, and after these preparation we compute the baryon mass shift 
of two-flavor baryons. Quantum state dependence of SU(2) instanton evaluated here is 
useful in Section [3j We also find that the mass shift of the baryon depends on its spin 
as well as the quantum number of its radial excitation, which was missing in [13]. As a 
consequence, the mass shift of the delta baryon, which has spin 3/2, is around 1.5 times 
larger than that of the proton and the neutron. This is a characteristic prediction of the 
Sakai-Sugimoto model, unlike conventional chiral soliton models. 

In Section [3J we compute a shift of the baryon mass due to the quark masses in the 
case of three flavors. A soliton solution of the massless three-flavor Sakai-Sugimoto model 
has been studied in [18], where a mass formula for baryons was obtained. The masses 
of SU(3) baryons are degenerate when quarks are massless. However, once the SU(3)- 
breaking quark masses are considered, there appear mass splittings of baryons. Using the 
SU(3) baryon configuration constructed by embedding the SU(2) instanton in a SU(3) 
gauge field, we obtain a formula of the mass shift of three-flavor baryons, which depends 
on the flavor charge, the radial excitation quantum number, and the spin of the baryons. 

Finally, in Section HJ we discuss the comparison of the mass splittings of baryons 
obtained numerically by our computation with those observed in experiments. In Table |3J 
we show the mass splittings of baryons with AY = 1 or AI = 1, while we show the 
mass splittings of baryons with A/3 = 1 in Table HI Agreements of the theoretical 
and experimental values in these two tables are qualitatively good but quantitatively 
inconclusive. The reason can be that the physical strange quark mass are no longer very 
small. To improve the matching, higher order corrections in the expansion in the strange 
quark mass are necessary. In Table El we show a mass difference of an octet baryon and 
its decuplet counterpart. This is presented for the study of the known problem concerning 
the overall magnitude of the baryon mass in the Sakai-Sugimoto model. We observe that, 
at the leading order in the quark masses, the inclusion of the quark masses does not 
improve the situation of this problem. 

2 A review of baryon mass shift 

For our computation of the mass shift of the baryons for three flavors in Section [3j here 
we give a brief summary of necessary ingredients, by reviewing [13J. 

Sakai-Sugimoto model 

The Sakai-Sugimoto model [U |2] is a five-dimensional gauge theory, in which Kaluza- 

problem in holographic QCD. In this paper, we simply ignore this problem. We expect that a certain 
decoupling limit while the QCD scale is fixed can be taken, once we know systematically the structure 
of the l/N c and 1/A corrections. 
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Klein(KK)-decomposed gauge fields are mesons, while baryons are provided by its solitons. 
The action is given by the following five-dimensional U(Nf) Yang-Mills- Chern- Simons 
theory in a curved background: 



S — Sym + Si 



>cs, 



S 



YM 



— K 



d 4 xdz tr 



'l + z 1 



i? 1 + 



1 + z 2 )7\ 



jlZ 



(2.1) 



and Scs is the 5-form Chern-Simons term. Here, fi,u = 0,1,2,3 are four- dimensional 
Lorentz indices, and z is the coordinate of the fifth dimension. The field strength is 
defined as = d^A v — dyA^ — i[A p , A v \ with the U (Nf) flavor gauge fields A^ and A z . 
They are decomposed as 



A 



A + A- 1 



y/2N f 



1 AT* 

A a T a + A- ' 



(2.2) 



/ 



where T„ 



Nf 



1) are generators of SU(Nf), and ljy. is the unit matrix of size 
Nf. We work in the normalization tr [T a T b ] = \8 a b- 

This model describes large N c massless QCD. There are two parameters in ( 12. ip : 
a mass scale Mkk for which we chose a unit Mkk = 1, and k = XN c /(21Qtt 3 ) with 
the 't Hooft coupling constant A. The parameter Mkk can be easily recovered by a 
dimensional analysis when we numerically evaluate the mass shift of the baryons. The 
two parameters are chosen as 

M KK = 949 [MeV], k = 0.00745 (2.3) 

to fit the experimental valued of the p meson mass m p ~ 776 MeV and the charged pion 
decay constant f w ~ 92.4 MeV. 



Quark mass in the model 

Quark mass can be introduced into the the Sakai-Sugimoto model through worldsheet 
instantons. Connecting the D8- and D8-branes corresponds to breaking the chiral sym- 
metry in the model, and the worldsheet instanton amplitude induces a quark mass term, 
which breaks the chiral symmetry at the Lagrangian level. This method has been first 
developed in [3] and jl], but these two papers differ in the regularization of the worldsheet 
instantons. In this paper, we follow [3]. 

The relation between the quark mass and the worldsheet instanton is as follows. Let 
us connect the D8-branes and the D8-branes by D6-branes. Then we can put a worldsheet 
instanton whose Euclidean worldsheet has the boundary defined by the color D4-branes, 
the D8-branes, the D6-branes and the D8-branes. The worldsheet instanton involves 
QlQr vertex of the left- and right-handed quarks living at the D8-D4 and the D8-D4 
intersections, respectively. This vertex is precisely a quark mass operator. 

3 If one rigorously treats the masslessness, numerical values in the chiral limit should be used for f n 
and Ttip . 
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In the gravity dual picture, where two of the corners on the D4-branes are smeared 
out by the background curved geometry, we can still put the worldsheet instanton in the 
same manner. The worldsheet instanton amplitude which includes boundary coupling to 
the Wilson line written in terms of the meson excitation A z on the D8-brane, is given by 



5S 



(Fx P tr 



M exp 



A z dz 



+ c.c. 



(2.4) 



which should correspond to the quark mass tern0. Here M is a quark mass matrix, and 
z = z m and z = —z m specify the location of the D6-brane on the D8- and the D8-branes, 
respectively. The constant subtraction — 1^ makes sure that (j2.4p vanishes when .4 = 0. 

This term f!2.4[) can be written in terms of a pion field U(x). The field A z relates to 
the pion field [H|2], 



P exp 



dzA z 



exp [2in(x)/f v ] = U(x) 



(2.5) 



where the decomposition of tt(x) is defined as n(x) = ir a (x)T a . Thanks to this relation, 
the term (12.41) can be rewritten as 



SS = J d 4 x 8L, 5L = ctr [M(U + U j - 21 



(2.6) 



Here, although Zyyi IS 9b finite value, we approximate the exponential in (12. 4 j) by ( 12. 5p . 



Meson mass 

The definition of U in ( 12. 5ft is a standard notation for the pion field in chiral pertur- 
bation theories, and it turns out that H 2 . 6 1) reproduces a well-known form of the quark 
mass term in the chiral perturbation theories [301 El]- A constant c appears as an overall 
factor in (12. 6p . but it is difficult to evaluate the constant c accurately because it comes 

4 This term breaks a part of the gauge symmetry that the original action (|2 . 1 j) has, in a way consistent 
with the explicit chiral symmetry breaking by the quark mass term. The reason is as follows. 

Let us consider for simplicity the case that all the Nf quarks share the same mass, that is, M is 
proportional to an identity matrix ljv/ ■ We have coincident D6-branes in this case. Because the parts 
z — z rn and z — —z m on the D8-brane are connected by the D6-branes, the gauge transformations 
there must be equal. Thus the gauge symmetry reduces to a restricted part (explicit breaking). Under 
this restricted gauge transformation, (|2.4p is invariant. Now, recall that the D8-brane gauge symmetry 
at z — > ±oo, i.e. at the boundary, corresponds to the chiral symmetry of the boundary theory [I] [2]. 
Therefore, if z m is close to z = oo, the restriction put on the gauge symmetry by the D6-branes is 
consistent with the explicit chiral symmetry breaking by the quark mass. The breaking pattern is such 
that the left and the right chiral symmetries are broken to their diagonal subgroup. 

In general, one can choose different masses for the quarks; M is not proportional to the identity matrix. 
The Nf D6-branes are not on top of each other, but are now at z = ±Zm where i — 1, 2, • • • , Nf. In 
this case, the gauge symmetry on the D8-brane is further restricted, and under this restricted gauge 
transformation (|2.4[) is invariant. The remaining gauge symmetry is consistent with the explicit breaking 
of the chiral symmetry (U(N f ) L x ?7(JV»r -> (Z7(l))^). 
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from sub-leading contributions of the worldsheet instantons in the curved space. As we 
shall see, using meson masses can be an alternative to the explicit use of the constant c 
and the quark masses. We show the cases Nf = 2 and Nf = 3 in turn. 

In the case that Nf = 2, ( 12. 6p can be parametrized by a pion mass m^. The quark 
mass matrix is given by M = diag(m u , m d ), where m u and m d are the up and down quark 
masses, respectively. They are related with the pion mass as 

ml = ^-(m u + m d ). (2.7) 
ft 

Note that charged and neutral pions have the same masses in this case, at the leading 
order. 

In the case that Nf = 3, we need three meson masses. The quark mass matrix is given 
by M = dia.g(m u ,md,m s ), where m s is the strange quark mass. Using the charged pion 
mass m^± , the charged kaon mass m K ± and the neutral kaon mass m K o , we have 

m l± = -^( m u + m d ), m 2 K ± = j^(m u + m s ), m 2 K0R0 = ^{m d + m s ). (2.8) 

J TV J TV J TV 

Since these are the simplest relations between quark and meson masses, we use these as 
input parameters. 

In the rest of this section, we work with Nf = 2. 



Bar yon in the model 

The baryon is identified as an instanton soliton localized in the four-dimensional x 
space (M = 1,2, 3, z) pQ. The instanton number of the Yang-Mills theory (12. ip corre- 
sponds to the baryon number. An explicit solution of the equations of motion of the 
action (12. ip has been obtained in [Tj3]. Here we collect the part we need. 

The part relevant to our computation is the A z component of the solution. The non- 
Abelian part is identical to the BPST instanton, which, in a singular gauge, is given 
by 

A z = - (x l - X% ee Abpst , (2.9) 

while the U(l) part is simply A z = 0. Here £ = \J \z — Z) 2 + \x — A| 2 , and Tj (i = 1, 2, 3) 
are the Pauli matrices. In (I2.9p . five of the eight moduli parameters of the instanton come 
in explicitly: p is the size of the instanton, and X M = (X 1 , X 2 , X 3 , Z) = (X,Z) is the 
location of the instanton. As well as these five, there are three SU(2) rotation moduli. 
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The classical value of the moduli is also obtained in [19] as^| 



^ = irv!' Zcl = ' (2 ' 10) 



We can put X 1 = without loss of generality. 



Mass shift of the baryons 

Now we compute the baryon mass shift, including the dependence on the quantum 
states of the baryon. At the leading order in the quark mass, a classical shift of the baryon 
mass is simply given by 

SM = - J d 3 x 5L[A% (2.11) 

where A cl is the classical solution (12. 9p . It receives corrections when we consider the 
quantum states of the baryon. 

We begin with the z integral in f 12 . 5 j) . As the A z configuration (12. 9 j) is proportional to 
a matrix xVj for any value of z, the path ordering in (12. 5 j) reduces to an Abelian problem 
[32]. Hence we obtain 

U = exp [«7(r)acVJ (2.12) 

with 

r = v^x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 , x l = — . (2.13) 

r 

Further, using (12. 12|) and the relation (12. 7p . we obtain a formula for the mass shift, 

8M = 4vr/2p3 m 2 x 1.104 . (2.14) 

Let us consider the quantum state dependence of the baryon. Each of the baryon state 
can be specified by the quantum numbers {/(= J), I3, n p , nz}, where / is the isospin, J 
is the spin, and n p and n z are quantum numbers associated with the moduli p and Z, 
respectively. However, the Z-dependence disappears because it comes in as z — Z in 
the baryon solution (12. 9p . and disappears in (j2.14|) . Dependence on the SU(2) rotation 
moduli also disappears because U + W is proportional to the unit matrix, and the SU{2) 
rotation U + W — >■ G(U + Z7')G' gives no effect. Therefore, there is no dependence on 
nz and J3, once / is specified. The mass of proton and neutron are equal here. 

5 The action (|2.1j) is at the leading order in 1/A expansion, while the size of the instanton (|2.10[) is 
©(A -1 / 2 ). In order to take sub-leading corrections into account, the authors of [12] have studied also 
the Dirac-Born- Infold (DBI) action, and they obtained the same classical value of the moduli as in (|2 . 1 0|) . 
However, the DBI corrections do not suffice all the 1/A corrections. In this paper, we simply assume (|2.1|) 
as our starting point. 



fir) 



71 



y/1 + p 2 /r 2 
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The mass shift, however, depends on the modulus p. We can show that the expectation 
value of p n is labeled by n p and the spin l/2(= J); we write this as (p n ) np ,i- To compute 
this, we use the eigenfunction for the quantization of p obtained in [19], 

R(p;n p ,l) = C npjl e- C »^ 2 p?'- 2 F(-n p ,/3f,c pP 2 ), (2.15) 

where A = 1 + + l) 2 + 4iV c 2 /5, c p = 16it 2 k/VE, and F(a,r,z) is the confluent 
hypergeometric function defined by F(a, 7; z) = Yl°k=o feufr w ^h (ct)fc = ct(a + l) •••(«+ 
fc — 1). Here C„ pj ; is the normalization factor, and we can normalize R(p ; n p , I) as 

POO 

/ dpp 3 i?(p;n p ,/) 2 = l, (2.16) 
for each n p and /. Let us consider the case that n p = 0. Then, Co,; is computed as 



Using this, we can compute the expectation value of p n as 



(p n )„ p = 0l i = I dp p 3+n R( P ; 0, if = ^; 2 ;;:; y - (2.18) 



r(A + n/2) 



' ' ' cTm 



It is straightforward to repeat this computation for other values of n p . For example, for 
n p = 1, we have 

,„„, _(,. n(n + 2) \ m + n/2) _( n(n + 2) \ 

(p K-H - [I + —^-) cTm) - [l + -^-) (P W (2-19) 

To obtain the baryon mass shift for a given quantum state of the baryon, we use the 
formula (" 12 . 1 8[) with n = 3 and rewrite k in terms of f w as k — nf 2 /4:, 

tp)»^-js[Mi) W) ■ (2 - 20) 

Using ( I2.20p . we obtain the mass shift for the baryon with n p = as 

^%M= 1.104 x 6 7 t m tT ] - P.") 

Although two parameters /„- and m p can be fit by (12.31) . only /„. is relevant for evaluating 
the baryon mass shift. If we use f n = 92.4 MeV and N c = 3, (I2.2ip becomes 

SM np= o, l=1 =41[GeV -i ]> (2>22) 



m 

6M. 



n —n I— 1 . 



where the first line corresponds to the mass shift of the proton and the neutron, which 
share the same mass shift in the two-flavor case, and the second line corresponds to that 
of the delta baryon. The worldsheet instantons in Sakai-Sugimoto model predict that the 
mass shift of the delta baryon is around 1.5 times larger than that of the nucleon and the 
proton, at the leading order in expansion in pion mass. The result (I2.22p is consistent 
with values obtained in chiral extrapolation of results of lattice QCD [TH [151 EH 02] • 

3 Three flavors 

Now we consider dependence on three-flavor quark masses by introducing the strange 
quark mass. When quarks are massless, the baryon solution in three-flavor Sakai-Sugimoto 
model has been studied in [15] . It was found there that the classical moduli are given in 
this case again by (12.101) (we can put X 1 = 0). As a correction to the massless case, we 
study the additional action given by (12. 6p with Nf = 3 to evaluate a shift of the baryon 
mass. 

Here we pay attention to how to construct the SU(3) baryon states, for we can compute 
other than that in a way similar to Section |2j The £77(3) baryon states, whose wave 
functions are written by the seven SU(3) rotation moduli of the twelve SU (3) instanton 
moduli, can be obtained in the following three steps; we first embed the SU(2) classical 
soliton solution into the SU(3) gauge field, next consider the SU(3) rotation moduli, 
and finally project this rotated-soliton configuration onto each baryon state. For this, 
instantons can allow the techniques developed for Skyrmions [27] (see also [33j [MJ [351 G23 



The SU (2) BPST instanton solution (I2.9P can be embedded in the ££7(3) gauge field 



This A z configuration is still proportional to a matrix xVj as in the case of two flavors. 
Therefore the path ordering in (12.51) can be evaluated in the same way as in Section [2[ 
As a result, we obtain 



where f(r) appearing here is the same as that appearing in (I2.13p . Here Uq represents 
the classical embedding of the SU(2) instanton. 

The classical embedding (13. 2p should be rotated so as to include the SU(3) baryon 
states. Note that although the actual moduli space is SU (3)/U(l), one can simply consider 
rotations in SU(3) space, with imposing a constraint on harmonic functions on ££7(3) [55] . 
Hence we can simply put the ££7(3)-rotated form 



EZ|). 



as 




(3.1) 




(3.2) 



U(x) = GU {x)G\ G e £77(3), 



(3.3) 
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in ( 12. 6p . and we obtain [27] 



5L 



Ac 

- — (1 - cos/(r)) 



(m u + m d + m s ) 



—^-{m d -m u )D^{G) 



2m e — m 



u-m d (8) 



where 



^?(G) = -tr(GtA a GA 6 ) 



(3.4) 



(3.5) 



is the Wigner's D function for the adjoint representation of 577(3), and A a (a = 1, . . . , 8) 
are the Gell-Mann matrices. Further, we use (I2.8P to write the mass shift as a function 
of meson masses. Then we obtain 



5M = Anfy x 1.104 x ± [(l - V^>$(G) - D§>(G) 



in 



K°,K° 



+ (l + V3Dg\G) - D$>{G)) m| ± + (l + 2/jg ) (G)) m* ± j . (3.6) 

This mass shift depends on the quantum states of 577(3) rotation and p moduli, while it 
does not depend on Z. 

Now we project the operators D^(G) and D^(G) onto each baryon state. The wave 
function of the SU (3) baryons is given by [27J [28] 

vMG) = v / dm77( D /; - ;2 /;-, /;:; ., 



(GO* = fcf-lj^'^JfG)*. (3.7) 



Here r is a label for the representation of SU(3); Y,I,I 3 are the quantum numbers of 
S'L r (3)fl a vor, where I is the isospin, 7 3 is its third component, and Y is the hypercharge; 
J, J3 are those of SU (2) spin C 5?7(3) r ight, where J is the spin and J3 is its third component. 
A constraint Y R = N c /3 = 1 (for N c = 3) is imposed on the hypercharge of 5 , f/(3) ri ght- In 
the last expression of ( 13 .7p . fi and v represent {/, I 3 , Y} and {J, — J 3 , 1} in abbreviated 
form, respectively. Using (13. 7p . the expectation value of the D functions can be evaluated 



as 



(D a f(G)) 



B 



dG *%(G)D%\G)* b (G) 



r r. 



5^ l a fi n ) I b v v 



r r. 



(3. 



The summation is taken over all occurrences of the representation r in the product of the 
representation 8 and the other r. A table of the SU(3) Clebsch-Gordan coefficients 
is given, for example, in [ID] . The results are shown in Tabled] |27] . 



6 This is expected to come from the Chern-Simons term of the Sakai-Sugimoto model. A problem on 
the derivation of the constraint has been reported [TB], and here we simply assume this constraint. 



9 



Table 1: (D™(G)/I 3 ) B and (D^(G)) B 
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ft 
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Using Table [TJ to take into account the SU(3) baryon state dependence, we obtain a 
formula for the mass shifts in the baryon spectra, at the leading order in the quark masses 
which are rewritten in the meson mass squared, 



4vr/y x 1.104 x - (a Q m 2 K 

ajc° + a Km 2 K ± + a w m 2 ^) , (3.9) 



with the coefficients do, clk and a n listed in Table EJ This result (13. 9ft is symmetric under 
the interchange of two flavors: u -f* d, u •f-)- s, and d •f-)- s. This symmetry is manifest in 
Table d 

Table 2: The values of a , ax and a n for each baryon state 
8 P N A E+ E° E- 2° S- 



3 4 _9_ 3 11 8 4 8 
5 5 10 5 10 5 5 5 

4 3 _9_ 8 11 3 8 4 

5 5 10 5 10 5 5 5 
8864 4 4 33 

55555555 

10 A++ A+ A A - E*+ E*° E*" S*° E*~ ft" 
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We evaluate the expectation value of p 3 in ( I3.9p . We can evaluate its expectation value 
by using the eigenfunction obtained in the case of three flavors [18J, 



R(p; n p , (p, q),l)= C npMl e" c ^ 2 / 2 pPM-W*F(-n fi , P Ml ; c pP 2 ), (3.10) 
where rj = 8 is the dimension of SU(3), C n jp, q )i is a normalization factor, and 



'(n-l) 2 2N 2 8 

Ml = 1 + \j KJ - T J - + -Yt + s^ p2 + q2 + ZtP + ti+PQ) ~ l ( l + 2 )- ( 3 ' n ) 
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Here the (p,q) = (1, 1), / = 1 and the (p, q) = (3,0), I = 3 cases correspond to the octet 
and the decuplet, respectively. The eigenfunction is normalized as J dpp^R 2 = 1. Starting 
with (13.1 Op . we obtain the expectation value of p n in the n p = 



rQ9( p, g)i + n/2 ) 

2 T(/3(p,q)l) 



(P )n p =0,(p )3 ),i = n/2 — • (3.12) 



We find that it is necessary to replace fii of the two-flavor case with 0( p , q )i of the three-flavor 
case. Computations for n p ^ cases can be done in similar manners. 
Thus we obtain the mass shift for the baryon state with n p = as 

AM llfUx 63/4 F (/W + 3/2) 

OM B ,n p =0,(p,q),l = 1-104 X 



1 

x- (a m^ + a K m 2 K ± + a,mj ± ) , (3.13) 
o 

where /3( p , g )z is defined in (13. lip , and the coefficients a , ax and a„- are listed in Table [3 
This is our main result for the mass shift at the leading order in the quark masses from 
the chiral limit. Using = 92.4 MeV and iV c = 3, we obtain 

5M Binp=0 ,oct = g (a m 2 K0 + a K m 2 K ± + a n m 2 n± ) x 7.9 [GeV -1 ], 

5M Sinp=0 ,dec = g (ao^l-o + a K m 2 K± + a^ml ± ) x 9.5 [GeV -1 ], (3.14) 

for octet (spin 1/2) and decuplet (spin 3/2) baryons, respectively. In Section HJ we discuss 
a comparison with experimental data. 

We can evaluate the p-dependence also for the n p — 1 excitation, which includes 
the Roper excitation. We have the following ratio of the n p — 1 and the n p = 
cases: (p 3 ) np =i,( P , q )i/ (p 3 )n p =o,{ P , q )i = {/3(p, q )i + 15/ 4 )/Ap>^- Hence the mass shift s for the 
n p = 1 excited baryon states are given by multiplying the right-hand side of (I3.13P by 
{/3(p tq )i + 15/4)//3( Pi9 )/. As a consequence, 7.9 and 9.5 in (13. 14j) are replaced with 12 and 
14, respectively. 

Lattice QCD simulation with a light strange quark has not been carried out yet. We 
hope that lattice computations of QCD in the future, preferably with large N c , may 
reproduce our analytic computations (I3.13P and (13 . 14h . 

A comment on higher order corrections is in order. Our result (I3.13P is at the leading 
order in the quark masses breaking the flavor symmetry, and also in the 1/N C expansion. 
There are higher order corrections to the action of the Sakai-Sugimoto model. The higher 
order corrections in the 1/N C expansion are from the higher order string loops, while 
those in the expansion in the quark masses can be obtained by the worldsheet instanton 
amplitudes with higher instanton numbers [3J. It is important to compute such higher 
order corrections, to get closer to real QCD in the holographic approach. 
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4 Discussions 



We compare the mass splittings of baryons numerically evaluated from (13. 14ft with the 
observed baryon spectra in PDG [41 J. Note that the mass shift ( 13.14ft is obtained in 
perturbation at the leading order in quark mass from the chiral limit, and one should not 
expect our theoretical spectra to agree well, particularly with observed values for baryons 
containing strangeness. The comparison here is a first step to include contributions from 
masses of three-flavor quarks. 

We use the following experimental values as input parameters of meson masses: 



140 [MeV], m K ± = 494 [MeV], = 498 [MeV]. (4.1) 



Together with (12. 3p . we use five input parameters, all of which come from the meson 
sector. 

Table 3: Mass splitting of baryons with AY = 1 or AI = 1 



8 AmA_7v 




Am s o_ A 


Am=o_ A 


Am=o_ s o 


theory [MeV] 2.4 x 10 2 


4.8 x 10 2 


2.4 x 10 2 


3.5 x 10 2 


1.2 x 10 2 


experiment [MeV] 1.8 x 10 2 


2.5 x 10 2 


77 


2.0 x 10 2 


1.2 x 10 2 




10 


Ams*o_A« 


Aras*o_£*o 


Amsi-_ s ,o 




theory [MeV] 


1.8 x 10 2 


1.8 x 10 2 


1.8 x 10 2 




experiment [MeV] 


1.5 x 10 2 


1.5 x 10 2 


1.4 x 10 2 





The mass splittings of baryons with AY = 1 or AI = 1 are shown in Table El where 
the mass difference of two baryons B 1 and B 2 is written as 

Am Bl _ B2 = m Bl - m B2 . (4.2) 

The entry Amgo_ A corresponds to AI = 1, while other entries are for baryons with 
AY = 1. Here we adopt neutral baryons except for the fl~ baryon in order to avoid the 
effects of the electromagnetic interactions, which is not considered here. 

This result, at the leading order in quark masses, qualitatively captures a tendency 
of the baryon mass splittings, while quantitatively it is unsatisfactory. This could be 
because theoretical values are obtained by a linear extrapolation to the physical strange 
quark mass (70 MeV to 130 MeV), which is no longer very small. To improve, we need 
higher order corrections, especially in the expansion in the strange quark mass. 

The mass splittings of baryons in an isospin multiplet and with A/3 = 1 are shown 
in Table HI We again find a qualitatively good but quantitatively inconclusive agreement 
in the comparison with experiments. One of the reasons is the same as in the previous 
paragraph; although the A/ 3 = 1 mass splittings are proportional to a small value m^— m u , 
which is translated into the difference between m^± and m K o go, the computation here 
is at the leading order in the strange quark mass. In addition to this, we expect that the 
electromagnetic interactions ignored here can be effective for the A/ 3 = 1 mass splittings. 
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Table 4: Mass splitting between baryons with A/3 = 1 



8 


Am N _ P 


Am s -_ s o 


Am s o_s+ Am=-_=o 


theory [MeV] 


2.1 


5.1 


5.1 8.2 


experiment [MeV] 


1.3 


4.8 ± 0.1 


3.3 ± 0.1 6.9 ± 0.3 



10 


A777./\baryons 




Ams*o_s*+ Am=.-_=.o 


theory [MeV] 


3.1 


3.1 


3.1 3.1 


experiment [MeV] 


(< 2) 


3.5 ± 1.5 


0.9 ± 1.4 3.2 ± 0.9 



Table 5: Mass difference of an octet baryon and its decuplet counterpart 





B* - B 


10-8 


A - N S*° - S° 


_ ^0 


5Mg* 

Arris* - 
ArriB 


- 8Mb [MeV] 

-b\ theory [MeV] 

*-s|ex P [MeV] 




3.7 x 10 2 
2.3 x 10 2 


4.6 x 10 2 1.6 x 10 2 
8.2 x 10 2 5.2 x 10 2 
2.9 x 10 2 1.9 x 10 2 


2.2 x 10 2 
5.9 x 10 2 
2.2 x 10 2 



Finally, we focus on octet-decuplet mass difference. A problem is known that the 
overall magnitude of the baryon spectra is not satisfactorily reproduced in the Sakai- 
Sugimoto model, once we use (12 .3p [TBI 119] - For instance, in the case of massless three 
flavors, the mass difference between octet and decuplet baryons was calculated in [TS] as 

M 10 - M 8 = 0.386208 x M KK ~ 3.7 x 10 2 [MeV] = M , (4.3) 

where Mrk = 949 MeV is used. However, the observed mass difference between the 
average value of octet baryons and that of decuplet ones is around 230 MeV. 

Once we include the quark mass, this mass difference is modified as shown in Table [5j 
where B and B* label octet and decuplet baryons, respectively. 5Mb* — 5Mb is the 
contribution from the nonzero quark masses, and the octet-decuplet difference is evaluated 
as 

ATO B *- B |theory = M + 5M B * - 8M B . (4.4) 

From Table El we observe that, at the leading order in quark masses, the inclusion of quark 
masses does not improve the situation of this baryon mass problem. We believe higher 
order corrections other than the quark masses are relevant for resolving this problem. 
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